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Abstract 

A representation of solutions of the wave equation with two spatial 
coordinates in terms of localized elementary ones is presented. Ele- 
mentary solutions are constructed from four solutions with the help 
of transformations of the affine Poincare group, i.e., with the help of 
translations, dilations in space and time and Lorentz transformations. 
The representation can be interpreted in terms of the initial-boundary 
value problem for the wave equation in a half-plane. It gives the solu- 
tion as an integral representation of two types of solutions: propagat- 
ing localized solutions running away from the boundary under differ- 
ent angles and packet-like surface waves running along the boundary 
and exponentially decreasing away from the boundary. Properties 
of elementary solutions are discussed. A numerical investigation of 
coefficients of the decomposition is carried out. An example of the 
field created by sources moving along a line with different speeds is 
considered, and the dependence of coefficients on speeds of sources is 
discussed. 
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1 Introduction 



The construction of exact representations of solutions of the wave equation 
in terms of elementary localized solutions is the subject of the present paper. 
Fourier analysis yields a representation of solutions in terms of plane waves. 
An analogue of Fourier analysis in our considerations is continuous wavelet 
analysis. We use here a special version of continuous wavelet analysis based 
on the transformations of the affine Poincare group [lj (Wavelets based on the 
Poincare group was discussed also in [2], [3j.) The representation obtained by 
us is associated with the initial-boundary value problem in the half-plane. It 
is given in terms of the affine Poincare wavelet transform of time-dependent 
boundary data on a line, which is the boundary of the half-plane. The wavelet 
transform is efficient in processing of functions of time and coordinate, which 
are results of observations and which have a multiscale structure pQ. The 
result of propagation of data is given by us in terms appropriate for data 
processing. 

The representations of solutions of the non-stationary wave equation as 
superpositions of elementary localized solutions have been presented for the 
first time in the work of Kaiser [4] . They were based on the analytic wavelet 
transform, and he used very special spherically symmetric elementary solu- 
tions. Representations of solutions in terms of elementary solutions from a 
wide class were proposed in [5], [6j, [7] and discussed in [8]. That represen- 
tation was based on continuous wavelet analysis. Elementary solutions are 
constructed from two chosen ones by means of a similitude group of trans- 
formations that contains shifts, scalings, and rotations. In present work, a 
solution of the non-stationary boundary problem for the wave equation has 
been decomposed in localized elementary solutions by means of the space- 
time version of wavelet analysis. Elementary solutions are constructed from 
four chosen ones by means of an affine Poincare group of transformations that 
contains shifts and scaling both in space and in time and the Lorentz trans- 
formations. Preliminary results were reported in [9j. The affine Poincare 
wavelet transform (APWT) is a coefficient in the decomposition of solutions. 
The region of integration is determined by the values of parameters under 
which the APWT is above a certain level. Numerical implementation and 
examples of calculation of APWT were discussed in our work [TO] . 

For the sake of simplicity we consider here only the case of two spacial 
dimensions. 

In section [2j we recall affine Poincare wavelet transform theory. In section 
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|3j we develop Poincare wavelet theory for solutions and construct representa- 
tions of solutions of the wave equation based on it. Then we give examples of 
mother solutions and discuss properties of elementary solutions obtained from 
them. The implementation and interpretation of Poincare wavelet transform 
are considered in section [5] with several examples. 

We suggest to determine mother wavelets by means of their Fourier inte- 
grals which has not explicit form in the coordinate domain. However we may 
take as mother solutions exact particle-like solutions found in [TT] , [12] , [13] , 
[T4] which are given by simple explicit formulas both in the coordinate and in 
the Fourier domains. The review of another exact localized solutions is done 
in [T5] . [16J. Solutions obtained by Lorentz transformations are discussed in 

m 

A representation obtained of solutions can find applications in the non- 
stationary diffraction theory, in the propagation and scattering of multiscale 
fields, for example, in seismology. We hope that our results may be useful 
also in astronomy to determine velocities of relativist ic objects or to deal 
with time-dependent problems of scattering. 

2 Basics of affine Poincare wavelet analysis 

First steps in the application of wavelet analysis methods are the definition 
of a class of functions % considered, a choice of a mother wavelet i/j e %, 
and the construction of a family of wavelets. The mother wavelet is a fixed 
function. A family of wavelets is functions obtained from the mother wavelet 
by means of transformations of a group. The set of wavelets should be dense 
in the chosen space %. Then it should be proved that any function / £ T~L 
can be represented as an integral superposition of wavelets. 

It is useful to consider together functions and their Fourier transforms. 
Below we introduce the notation 

--(t> 

and define the Minkovsky inner product as follows: (<r, x) — —ut + k x x. The 
method of wavelet analysis itself is not associated with the wave equation. 
We apply it to the boundary data function f(x) = /( c ^ x ) dependent only on 
one spatial coordinate x and time t.^The Fourier transform of any f(x) £ L 2 

1 Here c is a parameter of the dimension of speed. We use the notation, taking in 
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reads 

/V) = J d 2 X f(x)e-^. (2) 

R 2 

Lorentz transformations enables one to connect the coordinate x and the 
time t in a stationary coordinate system with the coordinate x f and the time 
t' in a moving one with speed v. It reads 

= MJ). ( 3 ) 

/ cosh (f) - sinh \ ^ _ ^ 
y -sinh^ cosh^ y ' an ^~ ~ c ' ' ' 

where is called the rapidity and i> is the speed of the moving frame. The 
frequency and the wave number in the moving frame and in the stationary 
frame given by vectors cr f and <r, respectively, are related as cr f = A_^cr. 

We are going to clarify the choice of the space H. We choose any ip(ct, x) e 
L 2 . Construct a family of wavelets with transformations of the affine Poincare 
group, i.e., with shifting, scaling and the Lorentz transformations: 

^(X) = ^(A-^) ,&=£), (5) 
li = {6,a,0}. (6) 

A wavelet ^(x) for a = 1, b = represents a wave packet in the stationary 
frame, such that it is the mother wavelet in the moving frame. The con- 
structed family is not dense in the whole L 2 . To show this, we analyze the 
Fourier transform of wavelets 

g (<t) = a^(aA^a)e-^ b) (7) 

and show that we cannot go beyond the region bounded by the lines uj = ±ck x 
under the action of Lorentz transformations. For example, a point from the 
domain k > 0,\k x \ < k cannot be mapped to any point in the domain 
k x > 0, |fc| < k x , by choosing a and (/). 

The explanation of this fact is below. Let the support of the Fourier 
transform of the mother wavelet in the moving frame be concentrated near 
a point ctq, i.e., ^(cr f ) = outside a neighborhood of a f . The support of 

mind further applications to the wave equation in two spatial dimensions where the wave 
number is k = |fe|, the frequency is u = ±ck and c is the speed of light. 
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this packet in the stationary frame lies in the neighborhood of cr , such that 
<r = a _1 A ( ^cr / . It is easy to check that the point (cjq/c, k x0 ) lies on the 
hyperbola (cjq/c) 2 — fc 2 = a~ 2 ((u) f /c) 2 — k x \). For any prescribed (c^o,/^) 
we can find a unique parameter a > 0. The parameter cj) can be found by the 
relation tanh(0 + / ) = c/c x0 /^o, where tanh(0 / ) = ck f x /uj f . As (j) increases, the 
point (c^o/c, k x0 ) tends to the asymptotes of the hyperbola uj = ±ck x . Let 
u)q > 0, c|A4 | < u)q. Then for any a and the point (cjq/c, k x0 ) lies in the 
same region: uj > 0, c|/c x0 | < ^o- 

Therefore we should decompose the space L 2 (IR 2 ) into a direct sum: 

4 

L 2 (IR 2 ) = © % j taking into account their Fourier transforms. The space 

%j comprises functions fj £ L 2 (]R 2 ), which have the support of their Fourier 
transform in Vj (see fig[l]) : 

Mx) = j^y 2 j d 2 <rf(v)e^. (8) 

The mother wavelet is any function ipj(x) £ Uj that satisfies an admissibility 
condition: 

C j= /^ l^^ oo. (9) 

7 |(w/c) 2 -fc2| 

The family of wavelets is constructed by formulas (|5]). 
The affine Poincare wavelet transform Fj (//) of the function /j is defined 
as follows pQ: 

^•(M) = y ^/^Xffi /x = {6,a,0}. (10) 

We denote functions from L 2 by Latin lowercase letters and their wavelet 
transforms by Latin uppercase letters. It is not necessary to extract fj(x) 
from f{x) to calculate Fj{ji) ) as follows from the Plancherel equality: 

The main fact that yields the decomposition of solutions is the Plancherel 
equality for the affine Poincare wavelet transform [I]: 

(fi,9i)u = ^J diAFjMGM, (12) 
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(a) Domains V\ and T>2 (b) Domains V3 and V4 

Figure 1: Domains Dj that are supports of the Fourier transform of functions 
from j — 1, 2, 3, 4. The set of points obtained by Lorentz transformations 
from the point with (j) = is a hyperbola. Arrows on hyperbolas show the 
direction of point movement from = — 00 to = +00. 
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where 

oo 

da 



dfi= I d(j)j — I d 2 b. (13) 
o 



Equality (12) is valid for any fj,Qj E %j if the admissibility condition ^ 
is fulfilled. 

The reconstruction formula following from the Plancherel equality reads 
fi(x) = ^-JdfiF j (j i )^ ll (x). (14) 
See the [A] for more detail. Any fix) G L 2 can be represented as fix) = 

tfjix)- 

3 Decomposition of solutions 

In the present section we give exact decompositions of solutions of the wave 
equation based on the formulas discussed above. The decomposition is as- 
sociated with the boundary problem for the wave equation in a half-plane 
y > 0. — oo < x < oo: 

Utt = C 2 (U XX +Uyy) (15) 



u(ct,x,y) \ y=Q = f(ct,x), (16) 
u(ct,x,y) | t < = 0. (17) 

where c is the speed, t is time, x and y are spatial variables. The function 



f(ct,x) in (16) is assumed to be square integrable: 

J d(ct) J dx\f(ct,x)\ 2 < oo. (18) 

R R 

The desired decomposition is expressed in terms of solutions running away 
from the boundary under different angles. It was first reported in [9]. 

We define four subspaces of solutions j — 1, . . . , 4 of the wave equa- 
tion using the Fourier transform of the boundary data f{cr). Any function 
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u j(Xi v) ^ L 2 belongs to S)j if its Fourier decomposition reads 

Uj(x,y) = J d 2 af{a)e^^e^V^y ) j = 1,2, (19) 

Uj(x,y) = J d 2 af{a)e^e-V^ 2 y ) j = 3,4, (20) 

where k = uu/c is the wave number for any f(cr) G L 2 . The functions 
Uj(X) y)) J ' = 1? • • • ? 4, are solutions of the wave equation. If a function f(cr) 
decreases at infinity not so slowly to produce a classical solution, the functions 
Uj(x^y) satisfy the wave equation in the distribution sense. The solutions 
Uj(x^y)ij = 1j2, are propagating, the solutions Uj(x->y)->j = 3,4, are expo- 
nentially decreasing and associated with details of the spatial behavior of /, 
which are smaller than the wavelength. Note that the reconstruction formula 

4 

for the Fourier transform yields Uj(x-> 0) = f(x)- 

3=1 

In what follows we will use the Fourier transforms of solutions with respect 
to t and x denoted by Uj(a; y) : 

uji^y) = /V)eV^^, J = l,2, (21) 
u^a-y) =f(<r)e-V^yj = 3,4. (22) 

The Fourier transforms of solutions with respect to x and y dependent on t 
are denoted by Uj(ct] k X) k y ): 

uj{ct; k x , k v ) = f((T)e-^,j = 1,2, <r = (^kj+¥ y , k x ) T . (23) 

The solutions Uj(ct, x, y), j = 3, 4, increase exponentially for negative y. The 
Fourier transform of them with respect to x and y do not exist. 

To develop a procedure of continuous wavelet analysis in the subspace 
we should choose a solution G We call it the mother solution. The 
notation of the mother solution begins with the Greek uppercase letter and 
that distinguishes it from the mother wavelet in L 2 , which is denoted by the 
Greek lowercase letter. The mother solution can be reduced to the mother 
wavelet if y — 0. Therefore iftj(x) = ^j(Xi 0) should satisfy the admissibility 
condition ([9]). 

Then we construct a family of solutions: 
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Each solution represents a wave packet considered in the stationary frame, 
which is a shifted and scaled mother solution in the frame moving in the x 
direction. These solutions are transformed to wavelets from the family ^ 
on the line y — 0, i.e., 0) = where ipj^x) are wavelets. 



The decomposition of solutions of (15) in terms of mother wavelets read 



4 

u(x,y) = £^-(x,y) 5 u Ax,y) = ^fd^F^^ix.y), (25) 



F J ( / x) = /rf 2 X /(x)^(x,0). (26) 

It is a solution because it is a decomposition of solutions. It is reduced to 
the reconstruction formula for the boundary data / if y — and therefore it 
satisfies the boundary conditions. Each of solutions moves or decreases away 
from the boundary. We do not take into account solutions that come from 
infinity and increase in the positive y direction to satisfy initial conditions 



(17) 



The contribution of each solution is determined by the coefficient Fj, 
which is the wavelet transform of the boundary data. 



4 Examples of mother solutions 

Mother solutions can be taken from a wide class. We give examples of mother 
solutions, which behave as particles. As a starting point we take a solution 
that has the Fourier transform localized near the point k x = 0, k y = x : 

*(ct ] k x ,k y ) = expl- ll ' y 2 ^je-^, (27) 

where uj = c\k\. If t — 0, then the solution in the coordinate domain reads 

^■^-dj^p-^-sr)- <28) 

If t 7^ 0, the solution cannot be found explicitly. But if t is small enough, then 
^i(c£, x, y) « ^i(0, x,y — ct). It can be obtained by expanding k = + 
up to linear terms: uj = ujQ + c(k y — x) + . . . and substituting it into the phase: 

k x x + k y y — cjt — k x x + k y (y — ct) + . . . . 
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On this way, for small time we obtain a packet-like solution which is filled with 
oscillations of a spatial frequency k and which has the Gaussian envelope. 
It moves along the y axis with speed c. 

For large t the solution should be found numerically. Exact solutions 
having similar local behavior were found in [12] and were named 'Gaussian 
wave packets'. Solutions and their Fourier transforms were given by explicit 
formulas and were studied analytically. An investigation of these solutions 
from the point of view of continuous wavelet analysis was performed in [14J. 

To be a mother solution in fii, this solution should contain Fourier com- 
ponents with uj > and k y > 0. The last condition is fulfilled only approxi- 
mately if the central frequency of the packet k is larger than the width of the 
support of the Fourier transform in k y direction l/cry, i.e. k<j\\ = y/p » 1. 



To obtain the mother solution, we assume that formula (27) is valid only 
for k y > and \I/(ci; kx, ky) = if k y < 0. To make the Fourier transform 
smooth, we introduce an additional factor exp(— ^-)^k y = \/k 2 — k 2 . The 
mother solution ^i(ct, x, y) is of the form 



^i(ct,x,y) 




o\{k y -xf a]ki 



(29) 



2 2 
exp(— l/ky) exp (i(k x x + k y y — kct)) , 

k = \k\. Because of the factor exp (— the solution satisfies the admis- 
sibility condition. We take a variable of integration k = <Jk^ + k% instead 
of k y and write the Fourier transform of (29) in the form of (21). Then we 
obtain 



— k 

Ky 



(j\{k y — nf 



rr 2 h 2 



1 

k lt 



(30) 



K = V^ 2 ~ • Fig. 2a shows this solution in successive time moments. It 
represents a wave packet moving away from the boundary y = 0. It has a 
finite energy, because it is localized in an exponential way. The same solution 
for other choice of parameters in the case x<j\\ = ^/p = 2 is shown in Fig. 3a 



The mother solution in io 3 is constructed in a similar way. We interchange 
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5 10 15 20 5 10 15 20 



(a) Solutions in the moving frame (b) Same solutions in the stationary frame, 

0=1 

Figure 2: Real part of mother wavelet with parameters k = 16, o\\ = \/2? 
g j_ \/2 in successive time moments: ct = 0, ct = 7.5, ct 15. 
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5 10 15 20 



(a) Solutions in the moving frame (b) The same solutions in the stationary 

frame, 0=1 

Figure 3: Real part of the mother wavelet with parameters x = 4, cry = 1, 
= 2 in successive time moments: = 0, ct = 7.5, ct = 15. 



11 



k and k xi put iy instead of y, and obtain 



* ,( ' ;,) = ^P 2 (31) 

exp • e'^^X 

The solution \^ 3 (c£,x,y) is a wave packet running along the boundary and 
decreasing exponentially in the positive y direction. The Fourier transform 
of the solution comprises components with wavelengths in the x direction, 
which are smaller than the wavelength in vacuum. Figj4a| demonstrates the 
solution in successive time moments. The wavelet analysis of solutions from 
f) 3 and might be useful for analyzing of surface waves, but a more detail 
study of the problem goes beyond the scope of this paper. 




(a) Solutions in the moving frame (b) The same solutions in the stationary 

frame, <f> = 1 

Figure 4: Real part of the mother wavelet ^ 3 with parameters x = 4, a\\ — 1, 
a± = 2 in successive time moments: ct = and ct = 7.5. 

The mother solutions in io 2 and S)^ can be obtained as follows: ^(ct, x, y) = 
Wi(-ct, x, y), * 4 (c£, x, y) = V s (ct, -x, y). 

Now we discuss analytically the properties of the family of solutions ob- 
tained by means of Lorentz transformations. The mother solution ^i(ct, x, y) 
represents a wave packet moving in the positive y direction and localized near 
the point x = 0, y = ct. In the Fourier domain, it is localized near the point 
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cr = (x, 0), cjq = ex. Now consider ^i(cf, x f ) y), where x',y,t' are coordi- 
nates in the frame that moves with speed v and rapidity (f). The coordinates 
in the moving frame are connected with the coordinates in the stationary 
frame by transformations Assuming that t' is small enough, we obtain 
^i(cf, x', y) « ^i(0, x', y — ct f ) = ^i(0, cosh0(x — c£ tanh 0) , y + sinh^x — 
cosh 0ci). The position of the center of the packet x ,?/o in the stationary 



frame is x = citanh^, y = ctj cosh0. The packet (29) found in the moving 
frame is numerically small outside an ellipse 

For the sake of simplicity we assume that it is a circle a = a\\ = a^: 

(x') 2 + (y - ct 7 ) 2 < 2a 2 . (33) 

An example of such a solution is shown in Fig. 2a, Taking into account J3) 
we find the domain in the stationary frame corresponding to the circle (33) 
in the moving frame. It is an ellipse: 

cosh 2 cj)(x - x ) 2 + ((y - ?/o) - sinh cj)(x - x )) 2 < 2a 2 . (34) 

This ellipse in the main axes reads 

Xix 2 + X 2 f < 2a 2 . (35) 

where 

Ai, 2 = e ±(/) cosh0. (36) 
The directions of the main axes are 

ei,2 = ( ± f" ) - (37) 

The new coordinates x, y are connected with the old coordinates 
y — y by the relation 

j/J \y-yoJ Vl e il \ e 2\J 

where the columns of the matrix U are vectors ej,j = 1,2. 
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Figure 5: The support of the wavelet in the Fourier domain for various (j) 



The axes of the ellipse in the stationary frame are not orthogonal to the 
direction of propagation and to wave fronts (see Figj2]). 

Similar considerations are valid in the Fourier domain. The support of 
the Fourier transform of the solution is distorted under the action of the 
Lorentz transform. Examples of calculations for various (j) are given in Figj5} 
a discussion of Fig j5] see below. 

5 Numerical calculations of coefficients in the 
wavelet decompositions 



A decomposition of solutions is given by a fourfold integral (25) . Here we 
show that calculations of such an integral can be efficient for some boundary 
data. For realistic signals the domain of integration may be sparse. The 
domain of integration is determined by the numerical support of the coeffi- 



cient in decomposition (25) , which is the wavelet transform of boundary 
data. The wavelet transform is non-negligible if the supports of a signal (the 
domain where the signal is non-zero) and of a wavelet from the family ^ 
have an intersection. The supports of Fourier transforms of a signal and a 
wavelet should also have an intersection. 

The support in the Fourier domain is governed by a scale and a rapidity. 
Figure Figj5] demonstrates transformations of the support of the wavelet in 
the Fourier domain for different values of the rapidity (j). The support for 
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(j) = is a circle. Upon the Lorentz transformations, the support takes the 
form of an ellipse and its center moves along the hyperbola. Here there is 
an analogy with the transformations of the support in a spatial domain. A 
change in the scaling parameter a will lead to a change in the hyperbola, 
where the centers of supports are located. 

Shifts in the spatial domain b rule oscillations of an integrand in (10) 
. An increase in oscillations results in a decrease in the wavelet transform 
and vice versa. If the support of a wavelet for some parameters matches 
the support of a signal and b is taken in such a way that oscillations are 
minimized, then the wavelet transform has a maximum. 

Now consider an example. We assume that a function / represents a 
field that is generated in the plane y = by six groups of monochromatic 
point sources, which move in the x direction with different speeds in plane 
y = —5000. In every group, the speeds are distributed with respect to the 
Gaussian law with a = 0.01c, a is a dispersion of the distribution, c is the 
speed in the wave equation, c = 1. The sources have different frequencies. 
The sources in the jth group, j — 1, . . . , 6, are characterized by the frequency 
ujj and the rapidity 0j, which corresponds to the mean speed. We choose the 
parameters as follows: uj\ — 1, 0i — 0.4; uj 2 = 1, 02 = 0.7; u)s = 1, 03 = 0.5; 
cj 4 = 0.9, 04 = 0.3; u 5 = 0.95, 5 = 0.5; u 6 = 0.95, 6 = 0.4. The boundary 



function / as a function of x and t is shown in Fig, 6a, We show that 




(a) The function f(ct,x) (b) The scale-rapidity diagram 5i(a, (j)) 

Figure 6: An example of the scale-rapidity diagram. 

analyzing the field in the plane y = by means of the wavelet transform, 
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we can determine the parameters of the sources. The sources are far from 
the observation plane, and we do not calculate the wavelet transforms in 
domains ©3,4, because the contribution of the components propagating along 
the surface are negligible. 

We show that the wavelet transform enables us to find the frequencies 
and speeds of sources. We are not interested in the determination of the 
position of sources; for this reason we introduce a function S(a, 0), which 
reads 

S 1 (a,<l>) = ^J d 2 b\FM\\ ^ = {6,a». (39) 



Such a function, corresponding to the field in Fig. 6a, as a function of the 



scaling a and the rapidity (j) is given in Figj6b| The scale-rapidity diagram 
shown is calculated with the following wavelet parameters: x = 4, a\\ = 
2^/55 w 15, a± = 8, the calculation mesh is from —128 to 128 with step 0.5 
in both x and ct. We see six maxima in this figure, which correspond to each 
group of sources. The position of the maximum a, (j) contains information 
about the frequency lj — 1/a and the average rapidity of sources. (We 
recall that the speed v is connected with (j) by the relation v/c = tanh</>.) 
The domain of integration in the decomposition formula must contain values 
of the parameters a, (j) such that S is not small. If we are interested only in 
the investigation of sources of certain frequencies and rapidities, we need to 
take the corresponding domain of integration. 

A more complicated situation where the sources are masked by noise are 
considered by us in [ID]. It is shown therein that the contribution of a noise 
can be eliminated, because it is located at small scales on the scale-rapidity 
diagram. The parameters of sources are found. 
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A The reconstruction formula 



We prove the Plancherel equality (12) for the affine Poincare wavelet trans- 



form (10) to make the paper self-contained. For the sake of brevity, we 



assume that / = g. With account of (13) the Plancherel equality reads 

oo 

Jd<t>J^Jd 2 b |F i (6,a,0)| 2 = C i J d 2 x\f 3 ( X )\\ (40) 

R R2 R 2 

where Cj is defined by relation (|9j). 

Applying the Plancherel equality (11) to the inner integral with respect to b 
and taking into account the definition of Fj, we obtain 



d 2 b \Fj(b,a, 



d 2 a 



f(a) ^(aA_^er) 



(41) 



Here we use the fact that Fj is a convolution of two functions fix) an d 

^jfi{— x)- Then it is used that the Fourier transform of ipj^—x) is ^j/X* 7 )? 
which is calculated by the formula ([7]). 

Substituting (41) in the left-hand side of (40) and interchanging the 



order of integration, we obtain 



oo 



d 2 b \Fj(b,a,, 



(2*) s 



OC 

/V)fy # / 



da 
a 



d 2 cr I (a) I do I — 0j(aA_^cr) 

(42) 

The next step is an analysis of the inner integral on the right-hand side of 



(42) . For the sake of definiteness, we assume that j = 1. We transform the 



variable of integration in the following way: 



aA_ 



k cosh 6 + k x sinh < 



k sinh 6 + k x cosh < 



' p cosh(0 + O )\ , / cosh /N 

K p sinh(0 + 0o ) / ^ \sinh c 

(43) 

Here we use the notation: k = pcosh0 o , k x = psinh0 o , p > 0. Instead of 
a and we take new variables (j)' — (j) + 0o and p' = ap, (j)' G (— oo, +oc), 
p' G (0, oo). Then the vector variable cr f is chosen instead of 0', //. It takes 
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values in V\. The inner integral for any fixed a is transformed as follows: 



Id 



, 2 f \^j{<r')\ = 

^ J t '2_j t /2 — 



c 



(fc')' 



(44) 
(45) 



Formula (12) has a generalization, in which /, g are any functions from 



L 2 and the wavelet transforms of the functions / and g denoted by Fj and 
Gj, respectively, are calculated with different mother wavelets, say ( and ^, 
respectively. In this case, formula (12) is valid, but Cj is calculated as 
follows: 



\Mcy-kt\ 



< oc. 



(46) 



The reconstruction formula (14) may be obtained from (12) formally if 



we choose as a function g any sequence of functions that tends to the Dirac 
S function. To pass to the limit under the sign of the integral in 
must impose additional restrictions on the functions / and g. 



(12) 



we 



B Ellipse transformations 



Assume that the packet (29) is negligible outside the ellipse, i.e. if a 2 (x') 2 + 
f3 2 (y — ct f ) 2 > 1 in the moving frame. In the stationary frame this ellipse is 
distorted as follows: 

a 2 cosh 2 (j){x — Xq) 2 + f3 2 ((y — y Q ) — sinh <j>{x — x )) 2 > 1 (47) 

The solution is located in the ellipse which reads in the main axes 

Xix 2 + X 2 f < 1. (48) 

where Aj, j = 1, 2 are eigenvalues of the matrix of the quadratic form in the 
left-hand side of the inequality (47) which are as follows 

Ai, 2 = ^ 2 + /32 2 >CQSh2(/) ± X - yj (a 2 + /? 2 ) 2 cosh 4 - Aa 2 f3 2 cosh 2 0) . (49) 
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The directions of main axes read 
ei, 2 = 



/3 2 - a 2 cosh 2 <p-/3 2 sinh' ± J (a 2 + f3 2 ) 2 cosh 4 - Aa 2 /3 2 cosh' 

2f3 2 sinh <p 



(50) 

The new coordinates x, ^ are connected with old coordinates x — x , y — yo 
by the relation 



a// \y-yo/ \i e ii i e 2 

where the columns of the matrix U are vectors ej,j = 1,2. 
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